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ARABIAN  SEA  MIXED  LAYER  RESPONSE  TO 
1994-1995  OPERATIONAL  WIND  PRODUCTS 


1.  INTRODUCTION 

The  Arabian  Sea  experiences  strong  surface  forcing  from  the  winds  of  the  regularly  alternating 
northeast  (November  to  February)  and  southwest  (June  to  September)  monsoons.  The  northeast 
(NE)  monsoon  is  characterized  by  winds  of  up  to  10  m  s-1  from  the  northeast,  parallel  to  the  coast  of 
Somalia.  The  stronger  southwest  (SW)  monsoon  has  wind  speeds  of  up  to  18  ms-1  in  the  opposite 
direction,  with  the  strongest  winds  concentrated  in  July  over  the  western  Arabian  Sea  in  a  feature 
called  the  Findlater  Jet.  The  winds  influence  sea  surface  temperature  (SST)  by  driving  turbulent 
mixing,  which  entrains  deeper,  colder  water;  by  forcing  coastal  and  open  ocean  upwelling;  and  by 
forcing  horizontal  advection,  which  brings  cool  water  upwelled  along  the  Arabian  and  Somali  coasts 
into  the  central  Arabian  Sea  and  warm  water  north  from  the  equator. 

Surface  heat  fluxes  also  play  a  central  role  in  forcing  the  Arabian  Sea.  The  consequences  of 
this  forcing  are  not  simple  to  interpret,  as  a  feedback  is  possible  between  the  SST  change  and  the 
surface  heat  fluxes.  The  lower  SST  that  results  from  mixed  layer  deepening  can  result  in  greater 
surface  heating  because  of  the  SST  dependence  of  the  longwave,  latent,  and  sensible  components 
of  the  net  surface  heat  flux.  This  would  oppose  the  tendency  of  the  mixed  layer  to  deepen. 

Efforts  to  better  understand  the  forced  response  of  the  upper  ocean  in  the  Arabian  Sea  have 
often  in  the  past  used  ocean  general  circulation  models  (OGCMs)  (Kindle  1991,  McCreary  and 
Kundu  1989,  McCreary  et  al.  1993).  However,  the  lack  of  in-situ  data  has  made  it  difficult 
to  assess  the  realism  achieved  by  such  models,  especially  when  forced  by  surface  wind  and  heat 
fluxes  that  themselves  had  large  uncertainties.  Still,  because  of  the  energetic  large  and  mesoscale 
circulation  in  the  Arabian  Sea,  a  local  field  study  may  not  collect  the  data  needed  to  determine  the 
impact  of  larger  scales  on  the  local  evolution  of  the  mixed  layer.  In  this  study,  we  use  both  new 
in-situ  data  and  an  OGCM  to  address  these  problems. 

As  part  of  a  recent  cooperative  research  effort,  year-long  time  series  have  been  obtained  from 
a  mooring  in  the  central  Arabian  Sea  (Weller  et  al.  1998).  The  mooring  provides  the  records  of 
local  forcing  and  local  response.  An  OGCM  (the  Naval  Research  Laboratory  Layered  Ocean  Model 
(NLOM))  is  used  to  supplement  the  mooring  data  and  provide  information  about  the  larger  scale 
variability  in  which  the  mooring  was  located.  However,  because  the  model  circulation  depends  on 
which  of  several  available  wind  fields  is  applied  as  forcing,  we  also  examine  the  sensitivity  of  the 
OGCM  predictions  of  mixed  layer  depth  (MLD)  and  SST  to  the  choice  of  forcing  fields. 

In  doing  this  study,  we  had  three  goals:  to  examine  the  relationship  between  wind  forcing  and 
mixed  layer  dynamics  in  the  Arabian  Sea,  to  contribute  to  the  interpretation  of  the  upper  ocean 
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variability  observed  at  the  mooring  by  using  the  NLOM,  and  to  assess  the  predictive  capability  of 
the  NLOM.  The  OGCM  is  fully  described  in  Section  2.  A  variety  of  forcing  fields  (described  in 
Section  3)  are  used  to  force  a  suite  of  model  simulations  to  explore  the  causal  relationship  between 
the  forcing  and  the  upper  ocean  dynamics.  OGCM  experiments  using  different  wind  products  are 
described  in  Section  4.  Sections  5  through  7  discuss  the  OGCM  simulations  of  the  SST,  MLD, 
and  currents  at  the  location  of  the  moored  array.  The  surface  heat  fluxes  from  the  model  runs  are 
compared  to  those  observed  at  the  moored  array  as  a  check  on  the  realism  of  the  model  in  Section  8. 
Section  9  discusses  the  heat  balances  for  the  mixed  layer  and  uppermost  dynamical  layer.  Finally, 
Section  10  presents  our  summary  and  conclusions. 

Four  appendix  are  included  in  this  report  to  provide  more  detailed  information  on  specific 
aspects  of  the  NLOM.  Appendix  A  describes  the  oxygen-based  diapycnal  mixing  scheme  used  in 
the  NLOM  for  the  global  mixing  correction  scale  factors.  Appendix  B  provides  a  useful  list  of 
the  many  symbols  used  in  the  equations  throughout  the  report.  Appendix  C  presents  the  detailed 
derivation  of  the  heat  equation  used  for  the  heat  balance  results  reported  here.  Finally,  Appendix  D 
provides  a  glossary  of  acronyms  used  in  this  report. 

2.  THE  MODEL  OCEAN 

The  model  used  for  this  study  is  the  reduced-gravity  version  of  the  NLOM  with  modifications 
to  embed  a  surface  mixed  layer.  The  model  without  the  mixed  layer  has  already  been  well-described 
in  Hurlburt  et  al.  (1996)  and  Metzger  and  Hurlburt  (1996),  and  been  used  to  examine  the  Arabian 
Sea  circulation  by  Bruce  et  al.  (1994),  Keen  et  al.  (1997),  and  Young  and  Kindle  (1994).  The 
reader  is  referred  to  these  for  the  details  on  the  circulation  component.  In  the  modified  model 
used  here,  the  surface  mixed  layer  is  embedded  within  the  upper  layer  of  the  circulation  model  in  a 
manner  similar  to  the  “fossil  layer”  approach  first  developed  by  McCreary  and  Kundu  (1989),  and 
later  augmented  by  them  for  the  Indian  Ocean  basin  (McCreary  et  al.  1993,  MKM  hereafter).  The 
reader  is  referred  to  these  references  for  the  physical  rationale  and  layer  dynamics  of  this  approach. 

2.1  Layer  Model 

The  NLOM  has  prognostic  equations  of  motion  for  transport  (Vfc),  density  (pfc),  and  layer 
thickness  ( hk )  as  given  in  Metzger  and  Hurlburt  (1996).  We  follow  their  notation  throughout  this 
manuscript1.  Four  significant  changes  were  made  to  the  OGCM  to  accommodate  the  presence  of 
the  mixed  layer  and  improve  upon  the  OGCM  simulation. 

The  first  change  is  the  vertical  mixing  velocity  for  the  first  interface  (wi)  now  includes  a 
contribution  due  to  turbulent  entrainment  (u>kt)  whenever  the  mixed  layer  has  merged  with  the 
first  layer. 

Jo,  ^  =  0 

*  (  0)^  +  —  (dfc  —  k  >  1 

0  for  h\  >  hm 

max(0,  ojkt)  for  hi  <  hm 


(1) 

(2) 


The  definition  of  u>kt  is  given  further  below  ( cf  Section  2.2). 

‘Physical  quantities  which  are  obtained  from  a  climatology  will  be  denoted  by  a  caret  (  )  throughout  this  report. 
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The  second  change  is  the  use  of  a  monthly  varying  climatology  ( p j£)  in  the  relaxation  term 
of  the  density  equation  ap(pk  -  Pk)  (cf  Eq.  (A2)  in  Metzger  and  Hurlburt  (1996)).  Traditionally, 
the  NLOM  has  relaxed  the  layer  densities  to  an  annual  mean  density  based  on  Levitus  (1982). 
However,  from  NLOM  trial  simulations  with  the  embedded  mixed  layer,  it  was  found  necessary 
to  relax  the  layer  densities  to  a  monthly  varying  climatology  in  order  to  get  the  correct  seasonal 
variation  in  the  layer  1  temperature.  By  this  we  mean  the  ap  varies  each  month,  with  the  monthly 
mean  p^  interpolated  in  time  between  the  midpoints  of  the  months.  Here,  the  relaxation  time  scale 
(cr"1)  is  equal  to  the  length  of  the  month.  The  monthly  mean  p^  is  determined  using  monthly 
mean  layer  thicknesses  ( h obtained  from  multiyear  NLOM  simulations  during  the  spin  up  phase 
(cf  Section  2.3).  The  layer  densities  are  obtained  from  the  monthly  temperature  and  salinity 
climatologies  of  Levitus  et  al.  (1994,  1994a)  using  an  equation  of  state  with  zero  pressure  (Millero 
and  Poisson  1981)  to  be  consistent  with  the  NLOM  assumption  of  an  incompressible  ocean.  The 
annual  Levitus  climatology  is  used  to  supplement  the  monthly  climatology  at  depths  beyond  1000  m. 


The  third  change,  the  thermal  expansion  coefficient  (a),  has  been  given  a  temperature  depen¬ 
dence  to  accomodate  global  scale  modeling. 


Po  oT 

—  — (a\  +  2a2T). 

Po 

This  is  derived  from  the  assumed  equation  of  state 

p(T,  S )  =  ao  —  a\  T  —  a2  T2  +  as(S  —  So)? 


(3) 


(4) 


which  is  fitted  using  the  global  annual  Levitus  climatologies  for  temperature  T  and  salinity  S. 
Table  1  gives  the  parameter  values.  The  values  (ao,ai,a2,  as,  So)  are  obtained  by  using  the  same 
equation  of  state  as  above  and  then  applying  a  least  squares  fit  to  all  densities  with  salinities 
between  30  and  40  psu  and  temperatures  between  10°  and  25  °C.  Earlier  NLOM  studies  (Metzger 
and  Hurlburt  1996)  assumed  a  linear  dependence  of  layer  density  on  layer  temperature  to  obtain 
the  heat  flux  forcing  term  in  the  density  equation.  A  temperature-dependent  thermal  expansion 
is  included  here  to  accommodate  global  scale  modeling,  where  significant  density  deviations  would 
occur  due  to  the  zonal  temperature  variations.  The  climatological  temperature  for  the  layer  (T&) 
is  chosen  as  the  most  appropriate  value  to  use  for  the  temperature  dependence. 


Table  1  —  Parameters  for  Fitted  Equation  of  State 


Parameter  Value 


Po 

1023.0  kgm  3 

ao 

28.27  at 

a\ 

0.0872  at  °C~l 

0.0044  at  °C~2 

as 

0.7674  at 

S0 

35.0  psu 
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In  the  NLOM  interlayer  mixing  scheme  (Hurlburt  et  al.  1996),  the  weighting  factor  Qk  is 
parameterized  to  represent  the  climatological  cross-isopycnal  detrainment  rate,  by  basing  it  on 
annual  oxygen  saturation  values  at  the  mean  interface  depths  (Levitus  and  Boyer  1994b).  Here  we 
calculate  Qk  using  a  variation  of  the  original  diapycnal  mixing  scheme  (Shriver  and  Hurlburt  1997) 
to  remove  anomalous  results  that  occured  with  the  original  formulation.  The  details  are  given  in 
Appendix  A. 

2.2  Mixed  Layer  Model 

The  surface  turbulent  boundary  layer  that  is  embedded  within  the  OGCM  is  assumed  to 
be  “well  mixed.”  The  temperature  Tm  is  therefore  defined  to  be  constant  throughout  a  layer  of 
thickness  hm,  where  the  latter  has  a  minimum  thickness  of  a  few  meters.  The  layer  thickness  hm 
is  defined  as  the  lower  bound  of  the  turbulent  boundary  layer  and  is  hence  the  MLD.  The  SST 
and  MLD  prognostic  equations  are  consistent  with  the  NLOM  dynamical  equations  for  layer  1,  and 
their  notation  is  given  in  Appendix  B. 


SST: 


&I}n 

dt 


+  vi  ■  VTm 


max(0,  u)m) 
hm 


{Tm  —  Te)  + 


Q 


Po^phn 


+  •  (hmVTm) 


(5) 


MLD: 


dhm 

dt 


+  V  •  {hmv  i)  —  Lom  . 


(6) 


Note  that  the  horizontal  velocity  within  the  mixed  layer  is  assumed  to  be  the  same  as  the  first 
dynamical  layer.  The  vertical  mixing  velocity  is  given  by 

(  0Jkt  for  hm  <  hi 

Um  ~  |  OJkt  +  £*>1  for  hm  >  hi  , 

where  ojkt  is  the  rate  of  mixed  layer  deepening  or  retreat  as  discussed  below.  Although  hm  >  hi  is 
allowed  in  the  calculation  of  com  and  u>i,  hm  <  h  is  enforced  for  numerical  stability  at  the  end  of 
each  time  step. 

The  temperature  Te  is  that  of  water  entrained  into  the  mixed  layer  from  the  thermocline 
below.  Its  value  must  be  obtained  by  parameterizing  the  thermocline  in  terms  of  model  variables 
(Ti,  Tm,  T2,  hi,  hm,  and  h2).  After  testing  several  parameterizations  (Gallacher  1994,  Gallacher 
and  Peggion  1999),  we  found  that  a  simple  linear  model  yielded  the  best  combination  of  accuracy 
and  robustness  for  providing  an  entrainment  temperature. 

Te  =  Ti  +  -^q— (T2-T1)  (8) 

hi  +  5/12 


The  temperature  Tk  of  the  kth  dynamic  layer  is  calculated  from  the  density  of  the  layer 
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Tk  =  Tk 


Pk  ~  Pk 
Oi{Tk)p0 


(9) 


Here  p^  is  the  kth  layer  density  climatology,  and  Tk  the  corresponding  temperature  climatology. 
The  departures  of  layer  density  (pk)  from  the  monthly  mean  climatology  {pi-  and  Tk)  are  assumed 
to  be  solely  due  to  temperature,  as  there  is  no  explicit  treatment  of  salinity  in  the  NLOM,  and  there 
is  no  surface  forcing  applied  due  to  evaporation  minus  precipitation  (E  -  P).  The  contribution 
from  salinity  can  be  incorporated  into  the  model  in  the  future  when  the  fluxes  of  fresh  water  are 
better  known. 


The  rate  of  mixed  layer  deepening  or  retreat  a ikt  is  determined  using  a  modified  version  of  the 
Kraus- Turner  mixed  layer  model  (Kraus  and  Turner  1967,  Niiler  and  Kraus  1977,  Gallacher  and 
Peggion  1999).  The  turbulent  kinetic  energy  (TKE)  budget  is  assumed  to  be  one  dimensional  (the 
vertical  fluxes  are  much  larger  than  the  horizontal  fluxes)  and  stationary  ( de/dt  =  0).  The  terms  in 
the  TKE  budget  are  parameterized  using  variables  integrated  through  the  mixed  layer  and  surface 
fluxes: 


/ 


de  3  nu  u  h™Tn,  u  \ 

—dz  -  m 3u*  -  —b*w*hm  -  —  b'w'(-hm) 

-(milt*  -  6msb*w*hm  +  m5fulhm ) 

=  0. 


(10) 


The  first  term  on  the  right  hand  side  is  the  shear  production  of  TKE  generated  by  the  surface  wind 
stress  (r),  ui  =  \r\/p0  is  the  friction  velocity,  and  e  is  the  TKE.  The  second  term  is  the  production 
or  damping  of  TKE  by  surface  buoyant  forcing: 


= 


gtt{Tm)Q 

PoCp 


(11) 


The  third  term  is  the  increase  in  potential  energy  in  the  mixed  layer  due  to  entrainment: 

dhm 


-b'w'(-hm)  =  a{Tm)gATn 


dt 


(12) 


To  be  consistent  with  the  equation  for  Tm .  we  should  have  A Tm  =  Tm-Te.  However,  it  was  found 
in  practice  that  the  mixed  layer  became  too  deep  during  the  intense  boreal  summer  cooling  in  the 
southern  Indian  Ocean.  In  order  to  obtain  mixed  layer  depths  consistent  with  climatology  (Rao  et 
al.  1989)  we  had  to  use  values  of  A Tm  =  4°C.  This  is  consistent  with  the  results  of  McCreary  et 
al.  (1993). 


The  last  terms  represent  the  turbulent  dissipation  of  TKE.  The  first  term  in  the  dissipation 
parameterization  is  a  fixed  fraction  of  the  shear  production.  The  second  term  is  a  fraction  of  the 
buoyant  production  of  TKE,  thus 


8  = 


1,  for  <  0 
0,  for  6*  re*  >  0. 


(13) 


The  third  term  was  not  included  in  the  original  Kraus- Turner  model.  It  incorporates  the  Coriolis 
time  scale  into  the  dissipation  parameterization, 
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f  =  max(|/|,/*),  (14) 

where  /*  =  2.54  x  10-5  s-1  is  the  Coriolis  parameter  at  5°  latitude.  This  is  particularly  important 
for  deep  mixed  layers  (Garwood  1977)  and  to  provide  the  correct  asymptotic  depth  in  the  neutral 
limit  (Blackadar  and  Tennekes  1968). 

The  model  constants  were  determined  from  one-dimensional  simulations  of  several  idealized 
cases  of  wind-stirring,  heating,  and  cooling.  The  mixed  layer  model  was  further  tested  with 
one-dimensional  simulations  of  the  seasonal  cycle  of  mixed  layer  deepening  and  retreat  at  Ocean 
Weather  Station  Papa  (Gallacher  and  Peggion  1999).  The  model  was  also  tested  in  an  idealize 
two-dimensional  study  of  air-sea  interaction  during  coastal  upwelling  (Gallacher  and  Peggion  1999). 
The  values  of  the  model  constants  are  given  in  Table  2  (see  Appendix  B  for  notation). 


Table  2  —  Turbulence  Model  Parameters 


Parameter 

Value 

cp 

3988  J  kg-1  °C_1 

A  t 

24  min 

ATm 

4  °C 

f* 

2.54  x  10-5  s-1 

9 

9.8  ms-2 

Ki 

20  m 

5.0  m_1 

m\ 

m3 

7.5  m-1 

m5 

5.0 

m6 

0.15 

n 

1 

When  the  production  of  TKE  exceeds  the  dissipation  and  damping  of  TKE,  the  mixed  layer 
deepens  until  the  increase  in  potential  energy  from  entraining  colder  water  equals  the  available 
TKE  (ATKE).  The  rate  of  mixed  layer  deepening 

dhm 
Ukt~  dt 

_  heq  (15) 

2A  t  ’ 

is  computed  by  calculating  the  equilibrium  depth  heq  at  which  this  balance  is  achieved  from  Eq.  (10) 
at  the  latest  time  step,  and  evaluating  its  change  relative  to  the  mixed  layer  depth  at  the  previous 
leapfrog  time  step  /i“.  The  equilibrium  depth  is  given  by 

heq  =  h~  +  ATKE/(agATrah- ),  (16) 

where  energy  is  assumed  to  be  converted  solely  by  entrainment  of  water  that  is  ATm  colder  than 
the  mixed  layer.  When  the  production  of  TKE  is  less  than  the  dissipation  and  damping  of  TKE, 
the  mixed  layer  retreats  to  an  equilibrium  depth  given  by 

m3U*3  -  miu*3 

eq  f  6*™*  -  5me,b*w*  +  rn5/u*2 


(17) 
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The  mixed  layer  TKE  budget  is  assumed  to  be  stationary,  i.e. ,  production,  dissipation,  and  mixing 
are  in  balance.  To  be  consistent  with  this  assumption,  the  mixed  layer  is  defined  to  instantaneously 
adjust  to  the  equilibrium  depth  heg  rather  than  as  a  relaxation  to  this  depth  with  a  characteristic 
time  scale.  Finally,  a  minimum  MLD  of  =  20  m  is  imposed  in  order  to  maintain  numerical 
stability. 

2.3  Model  Basin,  Boundary  Conditions,  and  Initial  Conditions 

The  model  domain  and  boundary  conditions  are  identical  to  those  used  in  earlier  Indian  Ocean 
studies  (Young  and  Kindle  1994,  Bruce  et  al.  1994,  Gallacher  and  Rochford  1995).  A  version  of 
the  NLOM  is  formulated  for  this  study  with  three  active  layers  within  the  upper  550  m.  The  fourth 
lowest  layer  (encompassing  the  volume  below  ~  550  m)  is  infinitely  deep  and  at  rest,  as  required 
in  a  reduced  gravity  formulation.  The  initial  undisturbed  layer  thicknesses  are  given  in  Table  3, 
along  with  the  complete  list  of  model  parameters  used  in  all  the  experiments  reported  here. 

Table  3  —  NLOM  Parameters  for  Indian  Ocean  Simulations 


Parameter 

Value 

Const. 

k=l 

CM 

II 

*- 

II 

00 

Ah 

300  m2  s_1 

Ho 

100  m 

kh 

5000  m2  s~x 

Hk 

115 

200 

235  m 

ht 

60 

75 

75  m 

K 

500 

600 

1000  m 

Uk 

0.01 

0.01 

0.01  ms-1 

ck 

0.0 

0.0 

0.0 

0.05 

0.05 

0.05  kg  m-3 

The  model  layer  density  and  global  detrainment  climatologies  are  initialized  by  vertically  inte¬ 
grating  the  Levitus  climatology  for  these  initial  layer  thicknesses.  The  initial  velocity  field  is  then 
determined  by  running  the  model,  with  no  mixed  layer,  from  a  state  of  rest  to  statistical  equilib¬ 
rium  of  the  domain-averaged  kinetic  and  potential  energy  for  each  layer.  The  model  is  forced  with 
the  monthly  mean  wind  stress  climatology  of  Hellerman  and  Rosenstein  (HR)  (1983)  and  no  heat 
fluxes  during  this  phase.  The  model  was  spun  up  for  10  years,  and  density  (pk)  and  temperature 
(Tfe)  climatologies  were  constructed  from  the  monthly  mean  layer  thicknesses  of  years  8  to  10,  and 
the  global  detrainment  climatology  was  constructed  from  the  corresponding  annual  mean  layer 
thicknesses.  This  10  year  spin  up  serves  as  the  initial  state  for  the  subsequent  spin  up  of  the  mixed 
layer  experiments  with  the  different  wind  products. 

3.  SURFACE  FORCING 

3.1  Surface  Wind  Stresses 

We  examine  the  sensitivity  of  the  mixed-layer  response  to  three  wind  stress  products  from  two 
operational  weather  forecast  models  as  well  as  wind  stresses  from  a  well-known  monthly  climatology. 
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By  comparing  the  responses  to  these  different  wind  products,  the  sensitivity  of  the  mixed  layer  to 
short-term  wind  events  in  the  operational  wind  stresses  can  be  determined. 

For  the  climatological  experiments  we  use  the  monthly  surface  wind  stress  climatology  of  Heller- 
man  and  Rosenstein  (1983),  and  the  wind  velocities  from  the  Comprehensive  Ocean- Atmosphere 
Data  Set  (COADS)  (da  Silva  et  al.  1994).  The  COADS  wind  velocities  are  used  in  the  heat  flux 
forcing  (cf  Section  3.2)  and  are  at  a  nominal  reference  height  of  20  m  above  sea  level.  For  the  first 
operational  wind  experiment  we  use  a  1990-96  data  set  of  12-hourly  wind  velocities  at  1000  mb  pres¬ 
sure  elevation  from  the  European  Centre  for  Medium-Range  Weather  Forecast  (ECMWF)  model 
(ECMWF  Research  Department  1988,  1991,  ECMWF  Technical  Attachment  1993).  For  the  re¬ 
maining  experiments  we  use  1990-96  data  sets  of  12-hourly  surface  wind  stresses  and  10  m  wind 
velocities  from  the  Navy  Operational  Global  Atmospheric  Prediction  System  (NOGAPS)  model  at 
Fleet  Numerical  Meteorology  and  Oceanography  Center  (FNMOC)  (Hogan  and  Rosmond  1991). 

For  the  ECMWF  product,  we  use  the  standard  bulk  formula  to  construct  wind  stresses  (r) 
from  the  1000  mb  winds  (v) 

T  =  PaCdMv,  (18) 

where  pa  is  the  density  of  air  and  Cd  is  a  constant  drag  coefficient.  To  obtain  a  wind  stress  with 
a  magnitude  consistent  with  the  HR  wind  stresses,  we  construct  a  pseudo  stress  |v|v  from  the 
monthly  mean  climatology  of  these  winds,  and  scale  their  Arabian  Sea  average  to  be  similar  to 
the  HR  wind  stresses  at  the  height  of  the  southwest  (SW)  monsoon.  No  consideration  is  given 
to  the  height  of  the  1000  mb  surface  in  calculating  the  wind  stress.  This  scaling  is  validated  by 
comparing  a  1994-1995  time  series  of  the  ECMWF  wind  stresses  against  those  computed  from  the 
Woods  Hole  Oceanographic  Institution  (WHOI)  mooring  meteorological  observations  (Weller  et  al. 
1998).  When  this  same  method  is  applied  to  the  FNMOC  surface  streses,  we  find  the  wind  stress 
is  underestimated  at  the  peak  strength  of  the  SW  monsoon.  The  FNMOC  surface  stresses  are 
therefore  multiplied  by  AT  =  1.35  to  bring  them  into  better  agreement. 

The  time  series  of  the  wind  speeds  from  the  various  wind  products  compare  well  with  the 
moored  buoy  data  over  most  of  the  1994-1995  sampling  period  (Fig.  1).  (Throughout  this 
manuscript,  12-hourly  averages  of  the  mooring  record  are  used  to  be  consistent  with  the  high¬ 
est  frequency  of  the  wind  products.)  A  5-day  moving  average  is  applied  to  help  distinguish  the 
differences  between  the  wind  products.  Wind  speed  differences  with  respect  to  the  buoy  are  on 
the  order  of  2-3 ms”1  outside  the  SW  monsoon,  and  differ  by  up  to  2-3 ms"1  during  the  SW  mon¬ 
soon.  The  ECMWF  1000  mb  winds  consistently  overestimate  the  buoy  wind  speed  and  the  COADS 
winds  consistently  underestimate  outside  the  SW  monsoon  period.  The  former  is  likely  because  the 
1000  mb  pressure  height  is  usually  closer  to  the  ocean  surface  than  10  m,  and  would  have  a  greater 
strength  according  to  a  log(.z/*o)  wind  profile.  The  reason  the  COADS  wind  speed  underestimates 
the  buoy  wind  speed  is  likely  because  they  are  constructed  from  ship  observations  taken  at  a  nom¬ 
inal  reference  height  of  20  m  (da  Silva  et  al.  1994)  and  would  therefore  be  weaker  than  the  10  m 
winds.  The  wind  scale  parameters  (Table  4)  yield  wind  stresses  that  also  compare  favorably  with 
the  buoy  observations  (Fig.  2).  Note  that  adjustment  of  the  paCd  values  allows  these  constructed 
wind  stresses  to  be  in  better  agreement  with  the  buoy  wind  stresses  than  the  operational  winds 
with  the  buoy  winds.  The  largest  differences  in  wind  stress  with  respect  to  the  buoy  again  occur 
during  the  SW  monsoon. 


A  v  (m  s' 


ML  Sensitivity  to  1994-1995  Wind  Products 


Fig.  1  —  a)  Time  series  of  wind  speeds  at  the  WHOI  mooring  location  (15.5°N, 
61.5°E)  from  four  wind  products:  the  COADS  climatology  (COADS),  the  ECMWF 
1000 mb  winds  (ECMWF),  the  FNMOC  10  m  winds  (FNMOC  10m),  and  the  WHOI 
buoy  Vector-Averaging  Wind  Recorder  (VAWR)  data  (BUOY).  A  five-day  moving 
average  was  applied  to  the  operational  and  buoy  winds  to  smooth  out  the  data.  The 
minimum  wind  speed  used  in  the  heat  flux  calculations  is  indicated  for  reference,  b) 
Differences  between  the  wind  speeds  of  the  wind  products  and  the  mooring. 
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Fig.  2  —  a)  Time  series  of  wind  stresses  at  the  WHOI  mooring  location  (61.5°E, 
15.5°N)  from  five  wind  products:  the  ECMWF  1000  mb  winds  (ECMWF),  the 
FNMOC  10  m  winds  (FNMOC  10m),  the  FNMOC  surface  wind  stresses  (FNMOC 
Surface),  the  HR  climatology  (HR),  and  the  wind  stresses  computed  from  the  WHOI 
buoy  meteorology  data  (BUOY).  A  five-day  moving  average  was  applied  to  smooth 
out  the  data,  b)  Differences  between  the  wind  stresses  of  the  wind  products  and  the 
mooring. 
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Table  4  —  NLOM  Experiments 


Expt. 

Wind  Stress 

Wind  Velocity 

Resolution 

K 

PaCd 

Clim. 

HR 

COADS 

2°  &  1° 

1.00 

ECMWF 

ECMWF  1000  mb 

ECMWF  1000  mb 

2.5° 

1.00 

1.5  x  10"3 

FNMOC 

FNMOC  Surface 

FNMOC  10  m 

1.25° 

1.35 

FNMOC  10m 

FNMOC  10  m 

FNMOC  10  m 

1.25° 

1.00 

2.52  x  10“3 

The  good  agreement  between  the  wind  products  and  observations  is  evident  from  the  high 
values  of  the  correlation  coefficients  r  for  both  the  wind  speeds  and  wind  stresses.  These  corre¬ 
lation  coefficients  are  determined  using  daily  average  values  of  the  wind  speeds  and  wind  stresses 
(Table  5)2.  To  obtain  a  measure  of  the  predictive  skill  of  the  different  wind  products,  we  have  also 
included  in  Table  5  the  nondimensional  skill  score  (Murphy  1995): 


SS(x,y)  =  1  - 


rmse(x,y)2 

°{y)2 


(19) 


where  rmse(x,y)  is  the  root-mean-square-error  between  x  and  y,  and  a(y)2  is  the  variance  of  the 
reference  value  y  (i.e. ,  mooring  observations).  All  wind  products  are  interpolated  to  the  same  grid 
resolution  for  the  simulations,  which  is  at  2Ax  of  the  model  grid.  A  detailed  comparison  of  the 
time  series  of  these  different  wind  products  may  be  found  in  Weller  et  al.  (1998). 

Table  5  —  Correlations  and  Skill  Scores  for  Wind  Products 


(Wind  Prod,  Buoy) 

Wind  Prod 

r(r,  t) 

SS(r,  r) 

r(v,v) 

SS(v,v) 

COADS 

0.83 

0.53 

ECMWF 

0.92 

0.72 

0.94 

0.38 

FNMOC  10m 

0.92 

0.77 

0.93 

0.85 

FNMOC  Surface 

0.91 

0.78 

HR 

0.83 

0.36 

3.2  Surface  Heat  Fluxes 

The  COADS  monthly  mean  atmospheric  climatology  is  used  to  calculate  the  surface  heat  flux 
forcing  (da  Silva  et  al.  1994).  The  forcing  is  calculated  following  the  same  approach  as  McCreary 
and  Kundu  (1989).  The  heat  flux  at  the  ocean’s  surface  is  considered  to  be  a  sum  of  four  parts: 
incoming  solar  shortwave  radiation,  outgoing  longwave  radiation,  sensible  heat  flux  (Qs),  and  latent 
heat  flux  (Ql).  We  elect  to  follow  the  MKM  and  do  not  include  model  SST  dependence  in  the 
outgoing  longwave  radiation.  Instead  we  construct  a  net  radiation  (QR)  forcing  from  the  addition 
of  the  shortwave  and  longwave  radiation  climatologies.  Only  the  sensible  and  latent  heat  fluxes 
have  a  dependence  on  the  model  SST  in  the  NLOM  net  heat  flux: 


2  All  correlation  coefficients  reported  in  this  manuscript  are  found  to  be  statistically  significant  upon  applying  the 
“no  correlation”  test  described  by  Chu  et  al.  (1998). 
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Q  —  Qr  +  Qs  +  Ql- 


(20) 


The  sensible  heat  flux  depends  only  on  the  air-sea  temperature  difference  and  the  surface  wind 
speed.  The  latent  heat  flux  depends  on  these  as  well  as  the  air-sea  difference  in  specific  humidity. 
The  latter  requires  calculating  the  saturated  specific  humidity  ( qs )  at  the  sea  surface  temperature 
(Tm).  The  heat  fluxes  are  given  by  the  expressions 

Qs  =  Qs  Qpa  Pa  ma x(vmin,  Va) 

*(Ta-  Tm) 

Ql  =  Cl  L pam&X-(VminiVa)  {Oa  ~  Qs) 

x@{qs  -  qa)  <£ 

Cl  =  Clo  +  b(Tm  Ta) 

_  ees 

Qs  ~  Pa  -  (1  -  t)es 

(  19.33Tm  \ 

es  =  6.1078  exp  ^  +  273--j, 

where  the  symbols  are  defined  in  Appendix  B  and  the  parameter  values  are  given  in  Table  6.  We 
adopt  the  oceanographic  sign  convention  of  positive  heat  fluxes  downwards.  Note  that  the  empirical 
coefficient  Cl  is  assigned  an  increasing  value  with  air-sea  temperature  difference  to  simulate  its 
observed  effect  on  atmospheric  stability,  and  that  it  is  limited  to  a  range  of  reasonable  values 
0.001  <  CL  <  0.002.  The  COADS  climatology  is  used  for  both  the  temperature  and  specific 
humidity  of  the  atmosphere  in  these  expressions. 

The  surface  wind  speed  va  is  constrained  to  have  a  minimum  value  of  vmrn .  This  is  to  account  for 
nonzero  latent  and  sensible  heat  flux  as  the  wind  speed  goes  to  zero.  A  value  of  vmin  =  4  ms  was 
found  to  yield  the  optimal  agreement.  This  was  determined  by  computing  the  net  surface,  latent, 
and  sensible  heat  fluxes  at  the  WHOI  mooring  location  using  the  COADS  wind  speed  and  SST,  and 
then  comparing  the  time  series  against  the  monthly  COADS  heat  fluxes  at  the  same  location.  This 
threshold  value  for  vmin  is  consistent  with  an  analysis  of  data  from  the  western  equatorial  Pacific 

Table  6  —  Surface  Heat  Flux  Parameters 


(21) 

(22) 

(23) 

(24) 

(25) 


Parameter  Value 


Cs 

Cpa 

Clo 

b 

e 

L 

Pa 

Pa 

vmin 


0.001 

1012  J  kg'1  “C"1 
0.0015 
0.00033  °C 
0.62197 

2.44  x  106  Jkg-1 

1013  mb 
1.1275  kg  m-3 
4  ms-1 
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Ocean,  where  derived  bulk  transfer  coefficients  were  found  to  increase  rapidly  with  a  decreasing 
wind  speed  below  4  ms-1  (Greenhut  and  Khaksa  1995). 

4.  WIND  EXPERIMENTS 

Four  model  simulations  were  performed  to  examine  the  mixed  layer  sensitivity  to  wind  products 
(Table  4).  The  first  simulation  is  forced  by  the  HR  and  COADS  wind  products  to  see  the  Arabian 
Sea  response  to  the  seasonal  reversals  in  monsoon  winds  (referred  to  henceforth  as  the  HR- COADS 
Expt.).  The  second  simulation  used  the  ECMWF  wind  product  (ECMWF  Expt.  henceforth).  The 
12- hourly  interannual  variability  of  this  wind  product  is  the  major  difference  with  respect  to  the 
climatological  wind  forcing.  The  third  experiment  uses  the  surface  stresses  (henceforth  FNMOC 
Expt.).  The  fourth  experiment  (henceforth  FNMOC  10m  Expt.)  uses  wind  stresses  constructed 
from  the  10m  winds  in  the  same  manner  as  the  ECMWF  wind  stresses  (Section  3.1),  with  an 
appropriate  value  for  the  drag  coefficient  (Table  4).  The  latter  two  FNMOC  experiments  are 
performed  to  explore  the  sensitivity  of  the  mixed  layer  to  the  effects  of  atmospheric  stability  on 
surface  stress,  as  the  FNMOC  surface  stresses  include  the  effects  of  dynamic  stability  within  the 
atmospheric  planetary  boundary  layer.  The  10  m  wind  velocities  are  used  in  the  heat  flux  forcing 
for  both  FNMOC  experiments. 

The  models  are  spun  up  from  the  initial  state  ( cf  Section  2.3)  with  heat  fluxes  obtained  with  an 
active  mixed  layer.  For  the  operational  wind  experiments,  the  model  is  first  forced  using  a  monthly 
mean  wind  climatology  constructed  from  years  1990-96  of  their  wind  stresses  and  no  heat  fluxes  to 
bring  the  model  to  equilibrium.  They  are  then  run  using  using  the  12-hourly  winds  from  January 
1990  to  January  1996,  starting  with  the  0Z  hour  wind  of  January  15,  1990.  This  choice  of  calendar 
day  coincides  with  the  end  of  the  climatological  wind  cycle,  which  terminates  at  the  monthly  mean 
for  January. 

5.  SEA  SURFACE  TEMPERATURE 

The  NLOM  reproduces  quite  well  the  seasonal  variability  in  the  SST  as  given  by  the  mooring 
record  and  the  climatology  (COADS)  when  forced  by  the  climatological  and  operational  wind 
products  (Figs.  3a-b).  The  largest  differences  between  the  model  and  buoy  SST  are  at  most  1.5  °C. 
The  SST  decrease  during  the  NE  and  SW  monsoons  (Table  7)  closely  follows  the  deepening  of  the 
mixed  layer  (Fig.  6),  indicating  the  importance  of  entrainment  for  bringing  deeper  colder  water 
into  the  mixed  layer,  and  thereby  cooling  the  surface  waters.  It  is  noteworthy  that  the  NLOM 
hindcasts  a  comparable  peak  SST  value  to  that  found  in  mid-May  of  1995  in  the  mooring  record. 


Table  7  —  Arabian  Sea  Monsoons 


Monsoon 

Start 

End 

1994  NE  Monsoon 

Nov.  1 

Feb.  15 

1995  Spring  Intermonsoon 

Feb.  16 

May  31 

1995  SW  Monsoon 

Jun.  1 

Sep.  15 

1995  Fall  Intermonsoon 

Sep.  16 

Oct.  15 
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Fig  3  _  Time  series  of  SST  from  the  four  wind  experiments:  a)  climatological  winds  (HR-COADS),  b) 
ECMWF  1000  mb  winds,  FNMOC  surface  stresses  and  10  m  winds  (FNMOC),  and  FNMOC  10  m  winds 
(FNMOC  10m).  The  SST  from  the  WHOI  mooring  (Buoy)  and  from  the  COADS  climatology  (COADS) 
are  included  in  a)  and  b)  for  comparison.  Panel  c)  contains  the  SST  differences  between  the  NLOM 
solutions  and  the  mooring. 
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Fig.  4  —  Power  spectral  density  (PSD)  of  the  mooring  a)  sea  surface  temperature  (SST),  b)  mixed  layer 
depth  (MLD),  and  c)  net  surface  heat  fluxes  (S.  Flux)  in  reciprocal  days  (d-1).  A  high-pass  Butterworth 
filter  has  been  applied  to  the  12-hourly  averaged  time  series  of  these  mooring  observations  to  remove  signals 
with  periods  greater  than  40  days.  The  40-day  cycle  at  /  =  0.025  tT1  is  an  artifact  of  the  signal  filtering, 
and  the  maximum  frequency  is  1  cT1  because  of  the  12-hourly  frequency.  The  PSD  has  been  scaled  by  the 
maximum  value  for  the  frequency  range  ( PSDmax ). 

Not  surprisingly,  the  12-hourly  operational  wind  products  are  more  successful  in  producing 
a  model  SST  that  follows  the  higher  frequency  variability  in  the  mooring  record.  This  is  most 
noticeable  in  the  double  SST  peak  of  April-May  1995,  and  in  the  earlier  rise  in  SST  during  August- 
October  1995.  The  peak  SST  value  in  mid-October  1994  is  also  well  reproduced,  albeit  a  few  days 
early.  Of  the  operational  wind  experiments,  the  FNMOC  yields  the  best  SST  agreement  with  the 
buoy  as  it  has  the  smallest  \ASST\  (Fig.  3c)  and  the  highest  SST  correlation  (Table  8).  We  note 
the  FNMOC  and  FNMOC  10  m  experiments  yielded  SST  solutions  that  are  quite  close  to  each 
other  during  the  mooring  observation  period. 

To  provide  a  measure  of  the  agreement  between  the  model  and  buoy  SST,  we  computed  cor¬ 
relations  between  the  time  series  (Table  8).  They  reveal  that  a  comparable  agreement  is  obtained 
between  the  model  and  buoy  SST  as  between  the  climatological  and  buoy  SST.  The  correlation 
between  NLOM  SST  and  buoy  SST  is  noted  to  increase  with  the  increase  in  value  of  skill  score 
between  the  wind  product  and  the  mooring  wind  record  (Table  5).  The  very  high  SST  correlations 
are  because  the  low-frequency  seasonal  cycle  dominates  the  time  series.  The  buoy-climatology 
correlation  is  the  strongest  because  the  higher  frequency  oscillations  in  the  SST  solutions  from 
the  operational  wind  experiments  are  out  of  phase  with  the  oscillations  in  the  buoy  SST.  The 
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Fig.  5  —  Coherency  spectrum  between  the  12-hourly  averaged  buoy  observations  and  time 
series  from  the  four  NLOM  experiments:  climatological  winds  (HR-COADS),  ECMWF 
1000  mb  winds,  FNMOC  surface  stresses  and  10m  winds  (FNMOC),  and  FNMOC  10  m 
winds  (FNMOC  10m).  The  coherency  spectrum  between  buoy  and  climatology  (Clim.)  is 
included  for  completeness.  The  coherence  spectra  are  shown  for  a)  SST  (Cssr),  b)  MLD 
(Cmld),  and  c)  net  surface  heat  fluxes  (C3ftx),  and  were  computed  using  a  fast  Fourier 
transform  of  length  32.  A  log  scale  is  used  to  accentuate  the  differences  at  low  coherence. 
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other  significant  frequencies  can  be  determined  from  the  power  spectral  density  (PSD)  of  the  buoy 
SST  (Stoica  and  Moses  1997).  After  a  10th-order  high-pass  Butterworth  filter  has  been  applied  to 
remove  signals  with  periods  greater  than  40  days,  the  PSD  determined  with  a  Hanning  window 
shows  that  the  daily  (/  =  Id-1)  and  10-day  (/  =  0.1  d_1)  variations  are  the  next  most  important 
contributions  (Fig.  4).  By  computing  a  coherency  spectrum  between  the  mooring  record  and  the 
other  time  series  (Fig.  5),  it  is  seen  the  FNMOC  experiment  yields  a  NLOM-buoy  SST  correlation 
that  is  slightly  better  than  the  corresponding  correlation  from  the  HR-COADS  experiment  for  time 
periods  of  1.25  to  3  days,  as  expected  for  a  twice-daily  surface  forcing. 

The  FNMOC  SST  solution  is  in  closest  agreement  with  the  buoy  SST  because  the  FNMOC 
surface  stresses  and  FNMOC  10  m  winds  follow  the  mooring  record  most  closely  (Figs.  1  and  2). 
The  surface  stresses  follow  the  relaxation  in  winds  in  the  mooring  record  quite  well  in  late  August, 
and  the  FNMOC  10  m  winds  have  a  deeper  and  narrower  minimum  in  late  July  compared  to 
the  ECMWF  winds.  Specific  time  periods  where  the  NLOM  SST  solutions  agree  well  (and  less 
well)  with  the  buoy  SST  can  also  be  directly  related  to  the  agreement  between  wind  products  and 
the  mooring  wind  record.  During  periods  when  the  wind  product  is  weaker  than  the  buoy  wind 
observations  (Fig.  lb),  the  model  net  surface  heat  flux  is  greater  because  the  reduction  in  wind- 
speed  dependent  latent  cooling  increases  the  model  SST  relative  to  the  buoy  SST.  The  reverse  is 
also  true  when  the  wind  product  is  stronger  than  the  buoy  wind.  Similarly,  whenever  the  wind 
stresses  are  weaker  than  the  buoy  observations  (Fig.  2b),  the  amount  of  cooling  due  to  entrainment 
of  cold  water  through  wind  stirring  is  decreased,  thereby  increasing  the  model  SST. 

It  is  noteworthy  that  although  the  ECWMF  wind  speeds  are  up  to  50%  stronger  during  the  SW 
monsoon  than  the  buoy  winds,  the  SW  monsoon  ECMWF  SST  differs  only  slightly  from  the  SST 
obtained  using  the  other  wind  products.  This  is  in  stark  contrast  to  the  ECMWF  SST  being  much 
lower  than  the  SST  from  the  other  simulations  during  the  NE  monsoon.  This  reveals  that  different 
physical  processes  are  responsible  for  the  decrease  in  SST  during  the  NE  and  SW  monsoons. 

We  also  note  that  the  more  pronounced  monsoon  relaxation  and  resurgence  in  the  last  week 
of  July  as  depicted  in  the  FNMOC  wind  products  relative  to  the  ECMWF  wind  products,  did 
not  produce  any  evidence  of  a  narrow  minimum  or  maximum  in  the  FNMOC  SST  compared  to 
the  ECMWF  SST.  Given  the  success  of  the  NLOM  in  capturing  the  SST  double  peak  in  the  1995 


Table  8  —  Time  series  correlations  (r)  of  the  monthly  climatologies 
(Clim.)  and  model  experiments  (Expt.)  with  respect  to  the  buoy 
observations  (Buoy). 


Expt. 

r 

(Expt., Buoy) 

r(Expt.,Clim.) 

r(SST,SST) 

r(MLD,MLD) 

r(Q,Q) 

r(SST,Ta) 

r(Q,Q) 

Clim. 

0.96 

0.63 

0.73 

0.98 

1.00 

HR-COADS 

0.89 

0.40 

0.64 

0.94 

0.80 

ECMWF 

0.90 

0.65 

0.56 

0.94 

0.70 

FNMOC 

0.93 

0.84 

0.61 

0.94 

0.71 

FNMOC  10  m 

0.92 

0.77 

0.65 

0.95 

0.76 
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spring  intermonsoon,  we  conclude  that  dynamical  processes  other  than  just  wind  forcing  must  be 
responsible  for  the  SST  minimum  that  appears  in  the  last  week  of  July  in  the  mooring  record.  One 
possibility  is  a  short  term  and  rapid  change  in  the  solar  radiation,  air  temperature,  and  specific 
humidity,  as  the  NLOM  simulations  have  been  forced  with  monthly  climatologies  of  these  quantities. 

6.  MIXED  LAYER  DEPTH 

The  buoy  MLD  used  for  comparison  against  the  model  MLD  prediction  is  chosen  as  the  depth 
at  which  the  temperature  has  decreased  by  0.1  °C  from  the  SST  (Rudnick  et  al.  1997).  For 
comparison  we  have  also  included  a  corresponding  climatological  MLD  which  we  have  inferred  from 
the  Levitus  monthly  temperature  climatology  (Levitus  et  al.  1994)  for  the  same  0.1  C  criteria.  This 
climatological  MLD  fails  to  properly  represent  the  seasonal  MLD  cycle  indicated  by  the  mooring 
observations.  It  indicates  an  actual  shallowing  of  the  MLD  at  the  peak  of  the  NE  monsoon,  and 
seriously  underestimates  the  depth  of  the  mixed  layer  during  both  monsoons.  We  ve  included  it 
to  merely  illustrate  the  difficulty  in  validating  OGCM  mixed  layer  depths  over  basin  scales  against 
MLD  climatologies. 

The  NLOM  is  less  successful  in  simulating  the  mooring  MLD  than  it  is  in  simulating  the 
mooring  SST  (Fig.  6),  having  differences  of  up  to  60  m  during  the  NE  and  SW  monsoons.  The 
MLD  correlations  (Table  8)  are  lower  because  of  the  20  m  minimum  depth  imposed  within  the 
model  to  maintain  stability.  This  prevents  the  NLOM  from  reproducing  the  very  shallow  buoy 
MLD  during  the  intermonsoon  periods,  and  the  10-day  oscillations  that  strongly  dominate  at  high 
frequencies  in  the  latter  (Fig.  4). 

The  correlation  for  the  HR-COADS  MLD  is  noticeably  poorer  because  the  mixed  layer  continues 
to  deepen  during  the  second  half  of  the  SW  monsoon  when  the  buoy  record  indicates  a  shallowing. 
This  is  confirmed  by  calculating  the  correlations  for  just  the  NE  and  SW  monsoon  periods.  For 
the  NE  monsoon,  a  significant  correlation  is  found  for  the  NLOM  MLD  with  respect  to  the  buoy 
MLD  (r  =  0.74),  while  a  very  poor  correlation  is  obtained  for  the  SW  monsoon  (r  =  0.05). 


The  operational  wind  products  are  more  successful  at  capturing  the  mixed  layer  deepening 
during  the  NE  and  SW  monsoons.  This  is  evidenced  by  their  higher  correlations  (Table  5)  due 
to  their  greater  coherence  at  periods  of  1.5  to  2  days  (Fig.  5).  All  operational  wind  experiments 
yield  a  deeper  MLD  during  the  NE  monsoon  than  the  HR-COADS  solution,  and  a  shallower  SW 
monsoon  MLD  that  more  closely  tracks  the  mooring  record  (smaller  \AMLD\  in  Fig.  6c).  They  are 
quite  distinct  in  their  NLOM  MLD  solutions  during  the  SW  monsoon,  indicating  the  sensitivity  of 
the  MLD  to  the  differences  in  the  wind  products.  The  FNMOC  solution  provides  the  best  MLD 
prediction,  in  large  part  because  it  captures  the  MLD  retreat  in  the  second  half  of  the  SW  monsoon 
as  indicated  in  the  mooring  record.  It  is  because  of  the  large  decrease  in  FNMOC  surface  wind 
stress  to  values  in  close  agreement  with  the  buoy  in  mid- August  that  the  FNMOC  MLD  captures 
this  retreat  event. 

The  deeper  FNMOC  mixed  layer  during  January  1995  relative  to  the  ECMWF  MLD  is  at 
odds  with  the  increased  surface  heating  one  would  expect  from  FNMOC  10  m  winds  that  are  up  to 
2  ms"1  weaker  than  the  ECMWF  1000 mb  winds  (Fig.  lb).  The  deeper  FNMOC  mixed  layer  is  due 
to  increased  net  surface  cooling  that  is  made  possible  by  a  northwestward  current  advecting  warm 
equatorial  surface  water  to  the  mooring  site  (c/  Section  9.1).  This  northwestward  current  is  induced 
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Fig.  6  —  Time  series  of  MLD  from  the  four  wind  experiments:  a)  climatological  winds  (HR-COADS), 
b)  ECMWF  1000  mb  winds,  FNMOC  surface  stresses  and  10  m  winds  (FNMOC),  and  FNMOC  10  m 
winds  (FNMOC  10m).  The  0.1  °CMLD  from  the  WHOI  mooring  (Buoy)  and  from  the  Levitus  climatology 
(Levitus)  are  included  in  a)  and  b)  for  comparison.  Panel  c)  contains  the  MLD  differences  between  the 
NLOM  solutions  and  the  mooring. 
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in  the  first  week  of  January  by  the  FNMOC  surface  winds,  as  opposed  to  a  northeastward  current 
as  obtained  in  the  last  2  weeks  of  January  with  the  ECMWF  1000  mb  winds.  This  northwestward 
current  brings  more  warm  water  to  the  mooring  location,  leading  to  a  warmer  model  SST  than  in  the 
ECMWF  wind  experiment.  As  a  consequence,  greater  surface  cooling  occurs  and  a  deeper  mixed 
layer  develops  in  our  model.  This  dynamical  response  is  due  to  model  limitations  in  representing 
the  surface  heat  fluxes  and  upper  ocean  turbulent  mixing.  The  air  temperature  is  prescribed  to 
remain  fixed  when  it  would  warm  up  in  the  prolonged  presence  of  the  warmer  SST,  thereby  reducing 
the  surface  cooling.  The  contribution  of  horizontal  advective  heat  fluxes  is  not  included  in  the  1-D 
turbulent  kinetic  energy  of  the  modified  Kraus- Turner  model.  The  inclusion  of  such  a  contribution 
would  increase  surface  buoyancy  and  thereby  inhibit  mixed  layer  deepening. 

Of  the  four  numerical  experiments,  the  FNMOC  surface  wind  stress  experiment  best  represents 
the  observed  MLD  during  the  second  half  of  the  SW  monsoon  (Fig.  6).  This  is  because  the  FNMOC 
surface  stress  best  depicts  the  major  relaxation  of  the  SW  monsoon  during  the  mid-July  to  mid- 
August  period.  The  stronger  FNMOC  10  m  and  ECMWF  1000  mb  wind  stresses  during  the  second 
half  of  the  SW  monsoon  are  responsible  for  the  deeper  mixed  layer  at  this  time. 

7.  CURRENTS 

The  mooring  record  shows  the  influence  of  strong  flows  associated  with  mesoscale  geostrophic 
features  because  of  the  strong  influence  of  eddy  activity  at  the  mooring  site  (Weller  1999).  As  a 
consequence,  the  NLOM  layer  1  currents  do  not  accurately  represent  the  corresponding  mooring 
layer  1  averages  (Fig.  7)  formed  from  the  time  series  of  mooring  subsurface  currents  and  layer  1 
thicknesses  (Fig.  8).  Accurately  predicting  currents  at  a  fixed  location  in  a  region  of  strong  eddy 
variability,  such  as  the  Arabian  Sea,  has  traditionally  been  a  difficult  challenge  for  all  OGCMs. 
While  the  model  is  able  to  represent  the  more  deterministic  boundary  currents,  an  accurate  depic¬ 
tion  of  the  mesoscale  variability  in  the  interior  of  the  Arabian  Sea  would  require  higher  horizontal 
resolution  (1/8°  or  better)  and  the  assimilation  of  satellite- inferred  SST  and  sea  surface  height 
(SSH)  data  into  the  model.  Studies  have  demonstrated  the  increased  model  skill  at  hindcasting 
mesoscale  eddies  and  meandering  frontal  structure  in  nondeterministic  regions  when  assimilating 
satellite-inferred  SSH  into  the  NLOM  (Smedstad  et  al.  1994,  Smedstad  et  al.  1997,  Hurlburt  et 
al.  1996).  The  improvements  are  smaller  in  equatorial  regions  where  the  ocean  has  a  strong  deter¬ 
ministic  response  to  wind  forcing.  NLOM  studies  have  also  shown  the  importance  of  having  high 
horizontal  resolution  and  upper  ocean-topographical  coupling  via  baroclinic  instability,  in  order 
to  properly  resolve  the  eddy  scales  (Hurlburt  et  al.  1996,  Hogan  and  Hurlburt  2000).  Including 
realistic  bottom  topography  via  the  finite-depth  formulation  of  the  NLOM  would  therefore  aid  in 
this  regard.  Implementing  these  features  is  beyond  the  scope  of  the  present  study. 

8.  SURFACE  HEAT  FLUXES 

The  NLOM  heat  fluxes  (Fig.  9)  agree  more  poorly  with  respect  to  the  buoy  than  the  climatolog¬ 
ical  heat  fluxes  (Table  5).  The  operational  wind  solutions  perform  more  poorly  than  the  buoy  and 
climatological  solutions  because  they  have  negligible  coherence  for  2  to  3  day  oscillations  (Fig.  5), 
even  though  they  perform  better  than  climatology  and  the  NLOM  HR-COADS  solution  for  5-day 
oscillations.  This  is  because  the  PSD  for  the  buoy  surface  heat  fluxes  occupies  a  larger  range  of 
higher  frequencies,  and  has  a  local  minimum  at  5-day  oscillations  (/  =  0.2  d  ,  Fig.  4).  The  lack 
of  better  coherence  with  increasing  frequency  for  the  NLOM  solutions  points  to  the  importance  of 
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Fig.  7  —  Time  series  of  the  NLOM  layer  1  currents  at  the  WHOI  mooring  location  (61.5°E,  15.5°N) 
for  the  four  wind  experiments:  climatological  winds  (HR-COADS),  ECMWF  1000  mb  winds  (ECMWF), 
FNMOC  surface  and  10m  winds  (FNMOC),  and  only  FNMOC  10m  winds  (FNMOC  10m).  The  lengths 
of  the  sticks  correspond  to  magnitude.  Northward  and  eastward  flow  are  in  the  upward  and  rightward 
directions,  respectively. 

other  factors  in  addition  to  wind  speed  for  predicting  surface  heat  fluxes,  e.g.  higher  frequency  air 
temperatures  and  net  surface  radiation. 

The  low  heat  fluxes  of  the  spring  intermonsoon  (Fig.  9)  coincide  with  the  mixed  layer  depth 
shallowing  to  the  imposed  minimum  value  (Fig.  6).  The  model  assumes  complete  absorption  of 
solar  radiation  heating  within  the  mixed  layer,  and  neglects  the  solar  heating  that  could  occur 
beneath  the  mixed  layer  when  light  penetrates  to  greater  depths  under  optically  clear  conditions. 
This  neglect  could  produce  excessive  SST  warming  when  the  mixed  layer  is  at  its  minimum  depth, 
thereby  producing  low  NLOM  heat  fluxes  because  of  excessive  latent  cooling.  Allowing  for  a  depth- 
dependent  solar  radiation  profile  is  expected  to  improve  the  agreement  between  the  NLOM  and 
buoy  surface  heat  fluxes.  The  excessive  surface  heating  by  the  model  during  the  SW  monsoon  is 
because  of  the  too-cold  SST  during  this  time  and  is  due  to  the  large  temperature  difference  imposed 
at  the  base  of  the  mixed  layer  ( cf  Section  5). 

A  remarkable  outcome  is  that  the  same  annual  cycle  occurs  for  the  surface  heat  fluxes  despite 
the  major  differences  in  the  wind  forcing  (Fig.  9c).  This  suggests  the  SST  is  relaxing  towards 
the  COADS  surface  air  temperature  in  a  manner  indicative  of  a  Haney-type  boundary  condition 
(Haney  1971).  To  test  the  conjecture,  correlations  were  computed  between  model  SST  and  COADS 
air  temperature,  and  between  the  model  and  climatological  net  heat  flux  (Table  5).  The  r(SST,  Ta) 
values  all  show  large  constant  correlations  (0.94-0.95)  that  are  comparable  to  the  correlation  be¬ 
tween  the  SST  and  air  temperature  from  COADS  (0.98).  The  r(Q.  Q )  correlations  of  the  model  net 
heat  fluxes  with  respect  to  COADS  are  appreciably  smaller  (0.70-0.80),  and  are  comparable  to  the 
correlation  between  the  COADS  and  buoy  net  heat  fluxes  (0.73).  If  SST  is  a  consequence  of  the  net 
heat  fluxes  applied  at  the  surface,  a  larger  variation  in  r(SST,Ta)  values  would  occur  that  follows 
the  variation  in  r(Q,  Q)  correlations.  The  absence  of  this  relationship  supports  our  conjecture  that 
our  heat  flux  forcing  is  relaxing  the  model  SST  towards  the  prescribed  air  temperature. 


This  relaxation  becomes  evident  when  the  heat  flux  term  in  the  SST  equation  is  expanded 
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Fig.  8  —  Time  series  of  a)  layer  1  temperature  and  b)  layer  1  thickness  from  the  four  wind  experiments: 
climatological  winds  (HR-COADS),  ECMWF  1000  mb  winds,  FNMOC  surface  stresses  and  10  m  winds 
(FNMOC),  and  FNMOC  10  m  winds  (FNMOC  10m). 
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Fig.  9  —  Time  series  of  net  surface  heat  fluxes  from  the  four  wind  experiments:  a)  climatological  winds 
(HR-COADS),  b)  ECMWF  1000  mb  winds,  FNMOC  surface  stresses  and  10  m  winds  (FNMOC),  and 
FNMOC  10  m  winds  (FNMOC  10m).  Net  surface  heat  fluxes  from  the  WHOI  mooring  (Buoy)  and  the 
COADS  climatology  (COADS)  are  included  for  comparison.  Panel  c)  contains  the  heat  flux  differences 
between  the  NLOM  solutions  and  the  mooring. 
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into  its  components  (Section  3.2).  A  linear  relaxation  of  the  model  SST  to  the  air  temperature  is 
introduced  via  the  sensible  heat  flux.  The  latent  heat  flux  introduces  a  nonlinear  relaxation  between 
the  model  SST  and  the  air  temperature  through  the  humidity  difference.  Because  there  is  no  change 
in  Ta  and  qa  in  response  to  the  thermal  state  of  the  ocean,  i.e. ,  no  air-sea  feedback  mechanism, 
the  NLOM  SST  relaxes  towards  an  equilibrium  state  that  is  offset  from  the  air  temperatuie.  As  a 
consequence,  the  net  surface  heat  fluxes  follow  a  repeatable  annual  cycle  that  deviates  only  when 
there  are  large  changes  in  the  ocean  model  MLD  and  SST.  This  is  consistent  with  the  findings 
of  other  authors  (Seager  et  al.  1988,  Seager  and  Blumenthal  1994)  who  argue  that  specifying  the 
air  temperature  and  humidity  in  the  heat  flux  boundary  conditions  of  an  OGCM  is  equivalent  to 
imposing  the  SST.  This  is  a  constraining  feature  of  this  surface  heat  flux  forcing.  However,  it  is 
the  best  heat  flux  formulation  that  can  be  implemented  using  monthly  mean  atmospheric  fields. 
Better  formulations  that  include  the  effects  of  atmospheric  stability  (Fairall  et  al.  1996,  Seager  et 
al.  1995)  require  forcing  the  NLOM  with  atmospheric  fields  that  vary  on  time  scales  of  24  hours 
or  less,  or  fully  coupling  the  NLOM  to  an  atmospheric  general  circulation  model. 

9.  LAYER  HEAT  BALANCES 

To  determine  the  relative  importance  of  the  different  heating  mechanisms  in  influencing  SST 
and  MLD,  we  perform  here  an  analysis  of  the  heat  balances  within  the  mixed  layer  and  the  up¬ 
permost  dynamical  layer  of  the  NLOM.  In  particular,  we  attempt  to  address  the  importance  of 
two  mechanisms  for  influencing  the  SST  during  the  1994  NE  and  1995  SW  monsoons.  The  first  is 
how  variations  in  surface  wind  stress  can  cool  SST  through  entrainment  of  deeper  colder  water  via 
wind  stirring.  The  second  is  how  changes  in  the  wind  driven  currents  may  influence  SST  via  heat 
advection.  We  include  a  heat  balance  analysis  of  the  first  layer  to  confirm  the  findings  of  the  heat 
balance  analysis  for  the  mixed  layer,  and  to  explain  changes  in  the  heat  balance  that  occur  when 
the  mixed  layer  and  layer  1  have  merged.  The  four  different  heat  balance  terms  due  to  surface 
heating,  advection,  interlayer  mixing,  and  diffusion  are  defined  in  Appendix  C. 

Although  the  model  does  not  represent  the  eddy-dominant  velocity  field  observed  at  the  buoy 
location,  we  use  the  model  circulation  to  examine  the  potential  role  of  the  basin-wide  background 
circulation  in  the  upper  ocean  heat  balances.  While  it  is  possible  to  compute  layer-average  horizon¬ 
tal  heat  advection  from  the  mooring  observations  to  compare  against  the  NLOM,  we  have  chosen 
to  omit  it  until  the  model  also  includes  the  assimilation  of  satellite  observations  to  better  depict 
the  effects  of  mesoscale  variability.  The  heat  fluxes  from  the  FNMOC  10  m  solution  are  not  shown 
because,  aside  from  the  horizontal  advection  contribution,  they  are  essentially  the  same  as  the 
FNMOC  solution. 

9.1  Mixed  Layer  Heat  Balances 


The  heat  fluxes  for  the  mixed  layer  (Fig.  10)  show  that  appreciable  horizontal  (H.  Adv.)  and 
vertical  (V.  Adv.)  advective  cooling  occurs  from  late  in  the  spring  intermonsoon  (early  May)  to 
near  the  end  of  the  SW  monsoon  (early  September),  in  counter  balance  to  the  heating  provided  by 
the  surface  heat  fluxes.  The  constraining  influence  of  global  detrainment  {Q,kQk)  on  mixed  layer 
deepening  is  responsible  for  the  decrease  in  vertical  advective  cooling  when  the  mixed  layer  and 
layer  1  have  merged  in  early  July.  Vertical  advective  cooling  ceases  completely  when  the  mixed  layer 
undergoes  retreat  in  late  August.  The  appearance  of  the  intense  Findlater  Jet  (Figs.  1  and  2)  is 
responsible  for  the  offshore  advection  of  water  along  the  Oman  coast  and  the  consequent  upwelling 
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of  cold  water.  It  is  the  arrival  of  this  advected  cold  water  at  the  mooring  (Fig.  11)  that  is  responsible 
for  the  change  in  early  July  from  small  heating  to  strong  cooling  in  the  horizontal  advective  heat 
flux.  There  is  little  heat  flux  from  diffusion  because  the  approximately  uniform  NLOM  SST  and 
MLD  produce  very  weak  horizontal  gradients. 

The  stronger  advective  heating  and  cooling  events  during  the  NE  and  SW  monsoons  in  the 
operational  wind  experiments  are  a  consequence  of  the  larger  changes  in  the  current  circulation 
induced  by  these  winds  (Fig.  7).  With  the  ECMWF  winds,  the  horizontal  advective  cooling  is  twice 
as  strong  in  the  1994  SW  monsoon  as  in  the  1995  SW  monsoon.  This  is  due  to  a  weakening  in 
the  1995  SW  monsoon  ECMWF  winds  relative  to  1994,  which  gives  rise  to  diminished  advective 
cooling  because  of  a  weaker  current  and  shallower  MLD.  The  ECMWF  horizontal  advective  cooling 
in  the  1995  SW  monsoon  is  greater  than  in  the  FNMOC  and  FNMOC  10  m  wind  experiments 
because  of  the  more  steady  eastward  current  produced  in  August  by  the  ECMWF  winds  (Fig.  7). 
All  operational  wind  products  display  advective  heating  during  the  1994  NE  monsoon,  with  the 
FNMOC  solution  being  the  strongest,  in  contrast  to  the  negligible  contribution  obtained  with  the 
HR-COADS  winds. 

The  average  and  integrated  mixed  layer  heat  balances  over  the  NE  and  SW  monsoons  (Table  9) 
reveal  that  the  decrease  in  SST  during  the  NE  monsoon  in  the  wind  experiments  is  due  to  both 
surface  cooling  and  vertical  advective  cooling  produced  by  the  combination  of  wind  stirring  and 
negative  buoyancy.  During  the  SW  monsoon,  vertical  advective  cooling  is  responsible  for  the 
decrease  in  SST,  with  the  significant  horizontal  advective  heat  flux  in  the  second  half  of  the  SW 
monsoon  making  a  secondary  contribution  of  up  to  as  much  as  one-third.  This  strong  cooling  from 
entrainment  is  offset  by  the  large  surface  heating.  The  monsoon  averages  for  the  net  surface  heat 
flux  from  the  model  are  consistent  with  mean  values  calculated  from  the  buoy  observations  (Weller 
et  al.  1998).  The  model  removes  more  heat  from  the  ocean  during  the  NE  monsoon  than  buoy 
measurements  indicate,  and  contributes  less  heat  to  the  ocean  over  the  SW  monsoon  than  recorded 
by  the  buoy.  Note  the  sum  of  the  model  heat  fluxes  remove  almost  twice  as  much  heat  energy  from 
the  ocean  at  the  mooring  location  over  the  SW  monsoon  than  over  the  NE  monsoon. 

9.2  Surface  Dynamical  Layer 

The  heat  fluxes  obtained  for  layer  1  for  the  various  experiments  (Fig.  12)  are  in  consistent 
agreement  with  those  obtained  for  the  mixed  layer  heat  fluxes  (Fig.  10).  Significant  heat  balance 
contributions  occur  for  all  terms  but  diffusion  (Fig.  12).  Diffusive  cooling  is  larger  here  than 
within  the  mixed  layer  (Fig.  10)  because  of  the  greater  spatial  variability  in  layer  1  thickness.  The 
heat  changes  due  to  advection  events  make  the  strongest  contributions  in  the  heat  balances.  For 
example,  for  the  climatological  solution  (Fig.  12a)  the  peak  advective  cooling  in  layer  1  is  twice 
as  large  as  the  peak  advective  cooling  in  the  mixed  layer.  This  cooling  occurs  because  of  the 
sharp  temperature  fronts  of  layer  1  created  when  cold  water  from  layer  2  is  entrained  into  layer  1. 
Clearly,  the  advective  heat  fluxes  are  sensitive  to  changes  in  the  basin  scale  circulation.  This  is 
evident  from  the  operational  wind  experiments,  where  significant  heating  and  cooling  advection 
events  for  layer  1  are  found  to  occur  (Fig.  12b-c)  for  which  there  is  no  mixed  layer  equivalent. 
These  strong  advective  events  correlate  with  the  changes  in  layer  1  temperature  and  thickness 
(Fig.  8),  and  are  attributed  to  geostrophic  flow.  For  example,  in  the  ECMWF  case  (Fig.  8b), 
the  layer  1  thickness  remains  approximately  constant  during  the  October-November  time  period 
as  compared  to  the  more  rapid  variation  in  the  climatological  wind  experiment  (Fig.  8a).  This 
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HR-COADS 


- H.  Adv. 

- H.  Diff. 

V.  Adv. 
- S.  Flux 


ECMWF 


- H.  Adv. 

- H.  Diff. 

V.  Adv. 
- S.  Flux 


FNMOC 


- H.  Adv. 

- H.  Diff. 

V.  Adv. 
- S.  Flux 


Fig.  10  —  Time  series  of  heat  fluxes  within  the  mixed  layer  for  years  1994-1995  from  the  four  wind  experiments: 
climatological  winds  (Clim.),  ECMWF  1000  mb  winds,  FNMOC  surface  stresses  and  10  m  winds  (FNMOC),  and 
FNMOC  10  m  winds  (FNMOC  10m).  The  layer  heat  fluxes  shown  are  the  contributions  due  to  horizontal  advection 
(H.  Adv.),  horizontal  diffusion  (H.  Diff.),  vertical  advection  (V.  Adv.),  and  net  surface  heat  fluxes  (S.  Flux). 
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Table  9  —  Shown  for  the  different  experiments  (Expt.)  are  the 
various  heat  fluxes  for  the  mixed  layer:  advection  (Adv),  diffusion 
(Diff),  layer  mixing  (LMix),  net  surface  heat  fluxes  (Sflx),  and  the 
net  heat  flux  from  their  sum  (Net).  The  first  row  (Buoy)  lists  the  net 
surface  heat  flux  values  from  Weller  et  al.  (1998)  for  comparison. 


NE  Monsoon  SW  Monsoon 


Expt. 

H.F. 

Ave.  H.F. 
(Wm“2) 

Buoy 

Sflx 

-19.7 

HR-COADS 

Adv 

0.21 

Diff 

-0.18 

LMix 

-16.07 

Sflx 

-29.97 

Tot 

-46.01 

ECMWF 

Adv 

1.53 

Diff 

-0.49 

LMix 

-25.75 

Sflx 

-24.79 

Tot 

-49.49 

FNMOC 

Adv 

17.69 

Diff 

0.81 

LMix 

-36.82 

Sflx 

-24.47 

Tot 

-42.77 

FNMOC  10  m 

Adv 

10.92 

Diff 

1.88 

LMix 

-34.69 

Sflx 

-20.12 

Tot 

-42.01 

Int.  H.F.  Ave.  H.F.  Int.  H.F. 
(MJm"2)  (Wm-2)  (MJm-2) 


-182.12 

89.5 

827.41 

1.89 

-41.28 

-391.57 

-1.62 

0.01 

0.09 

-143.96 

-118.83 

-1127.29 

-268.56 

69.60 

660.25 

-412.24 

-90.50 

-858.52 

13.72 

-13.72 

-130.13 

-4.36 

-0.49 

-4.64 

-230.68 

-115.52 

-1095.91 

-222.07 

56.15 

532.67 

-443.39 

-73.58 

-698.01 

158.54 

-7.36 

-69.85 

7.30 

0.80 

7.56 

-329.85 

-124.18 

-1178.04 

-219.21 

66.80 

633.70 

-383.23 

-63.95 

-606.63 

97.85 

-40.68 

-385.89 

16.88 

0.65 

6.17 

-310.82 

-117.36 

-1113.35 

-180.31 

66.63 

632.11 

-376.40 

-90.75 

-860.96 
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Fig.  11  —  The  SST  and  surface  currents  predicted  for  the  Arabian  Sea  basin  by  the  NLOM  wind  climatology 
experiment  (HR-COADS)  for  Julian  day  191  (July  10).  The  location  of  the  WHOI  mooring  at  (61.5°E, 
15.5°N)  is  indicated  by  the  pink  cross.  Clearly  evident  arc  the  colder  surface  waters  along  Oman  and  south 
of  Socotra  produced  by  coastal  upwelling  of  deeper  water.  The  advection  of  the  Oman  coastal  water  to 
the  mooring  is  partially  responsible  for  cooling  the  mixed  layer  at  this  location  in  mid-July. 

produces  a  much  weaker  geostrophic  flow  and  hence  a  considerably  diminished  advective  heat  flux. 
Given  that  the  transport  in  the  NLOM  is  directly  proportional  to  the  gradient  in  layer  temperature 
and  thickness,  it  is  not  surprising  that  advective  heat  fluxes  correlate  with  the  layer  1  temperatures 
and  thicknesses. 


As  expected,  entrainment  cooling  of  the  mixed  layer  only  occurs  during  the  brief  time  the 
mixed  layer  has  merged  with  layer  1.  It  contributes  much  less  towards  the  cooling  of  the  layer 
as  advection.  This  is  in  contrast  to  the  mixed  layer  case  where  entrainment  cooling  dominates 
throughout  the  period  of  deepening  due  to  wind  stirring. 

Advection  and  entrainment  are  responsible  for  cooling  layer  1  during  the  SW  monsoon.  The 
heat  input  over  the  SW  monsoon  from  the  surface  fluxes  is  approximately  balanced  by  the  heat 
removed  by  these  processes.  The  cooling  of  layer  1  during  the  NE  monsoon  occurs  because  of 
surface  cooling.  The  surface  heat  flux  is  responsible  for  the  net  heating  at  the  mooring  location 
over  the  course  of  the  year. 

10.  SUMMARY  AND  CONCLUSION 

In  this  report  we  investigate  the  sensitivity  of  mixed  layer  properties  in  the  Arabian  Sea  to  winds 
by  using  an  OGCM  forced  with  different  wind  products  to  interpret  a  one-year  record  of  mooring 
observations.  This  was  done  to  achieve  a  three-fold  goal.  First,  to  determine  the  sensitivity  of 
the  mixed  layer  in  the  Arabian  Sea  to  variations  in  the  winds.  Second,  to  identify  the  origin  of 
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mesoscale  features  and  help  interpret  how  they  are  responsible  for  features  observed  in  the  mooring 
record.  Third,  to  determine  which  wind  product  works  best  as  forcing  for  an  OGCM  when  seeking 
to  predict  mixed  layer  properties.  These  goals  were  achieved  by  carrying  out  NLOM  simulations 
using  wind  products  from  monthly  mean  climatologies  (HR  and  COADS)  and  operational  weather 
forecast  models  (ECMWF  and  FNMOC). 

We  found  that  both  SST  and  MLD  prediction  in  an  OGCM  are  sensitive  to  variations  in  the 
wind  stress  magnitude,  and  to  a  lesser  extent  the  wind  speed  dependence  in  the  surface  heat  fluxes. 
The  ability  of  an  OGCM,  in  this  case  the  NLOM,  to  reproduce  the  mooring  record  of  SST  and 
MLD  is  directly  related  to  how  closely  the  wind  product  used  to  construct  the  surface  fluxes  for 
the  model  agreed  with  the  observed  wind  stress  and  wind  speed. 

Despite  the  failure  of  the  NLOM  to  reproduce  the  current  variability  at  the  mooring  location, 
this  study  uncovered  a  significant  result.  The  large  surface  cooling  observed  in  the  mooring  record 
during  the  NE  monsoon  is  made  possible  in  part  by  a  northwestward  current  horizontally  advecting 
warm  water  to  the  mooring  location.  Were  it  not  for  this  advection,  the  SST  would  decrease  more 
than  is  observed  and  thereby  diminish  the  amount  of  surface  cooling  as  directly  seen  in  the  ECMWF 
wind  experiment. 

Of  the  wind  products  considered,  the  FNMOC  surface  wind  stresses  and  FNMOC  10  m  wind 
speeds  had  the  strongest  correlations  with  the  mooring  wind  record,  and  therefore  yielded  the  best 
agreement  between  the  NLOM  predictions  and  mooring  record  for  SST  and  MLD.  This  is  because 
the  surface  stresses  included  the  effects  of  atmospheric  dynamic  stability,  and  the  10  m  winds  were 
at  a  constant  reference  height.  The  ECMWF  1000  mb  wind  stresses  performed  more  poorly  because 
these  effects  were  ignored  (a  constant  drag  coefficient  was  used)  and  because  no  consideration  was 
given  to  the  height  varying  1000  mb  surface.  The  ECMWF  10  m  winds  were  not  used  in  this  study 
to  form  the  wind  stresses  because  in  our  earlier  experience  we  found  them  to  have  serious  bias 
problems.  Our  success  with  the  FNMOC  wind  products  indicates  the  importance  of  including 
atmospheric  dynamic  stability  when  constructing  surface  fluxes  for  use  with  an  OGCM  having  an 
embedded  mixed  layer. 

While  this  Indian  Ocean  NLOM  has  shown  good  skill  at  reproducing  the  mooring  observations 
for  SST  and  MLD,  it  is  clearly  deficient  at  predicting  the  currents  and  needs  improvement  in  the 
surface  heat  fluxes.  As  we  discuss  in  Section  7,  to  improve  upon  the  basin  scale  circulation  it  will 
be  necessary  to  increase  the  horizontal  resolution  of  the  model  to  at  least  1/8  so  that  the  NLOM 
will  be  sufficiently  eddy-resolving.  Moreover,  the  finite-depth  version  will  need  to  be  implemented 
to  take  full  account  of  the  coupling  between  the  wind-driven  surface  currents  and  the  topography- 
driven  abyssal  currents,  as  revealed  from  other  NLOM  studies.  Reformulating  the  model  to  allow 
for  a  “floating”  mixed  layer  that  can  penetrate  beneath  the  first  layer  interface  would  also  be  a 
major  improvement.  This  reformulation  would  reduce  the  increase  in  geostrophic  flow  that  occurs 
during  the  SW  monsoon,  when  the  first  layer  interface  is  deepened  by  the  mixed  layer,  and  thereby 
improve  the  layer  1  currents.  This  would  also  remove  artifacts  that  appear  in  the  model  SST  and 
MLD  as  a  consequence  of  constraining  the  mixed  layer  to  within  the  first  layer.  A  “floating”  mixed 
layer  has  already  been  implemented  in  global  version  of  the  NLOM  and  is  presently  being  evaluated. 

Immediate  improvement  in  the  surface  heat  fluxes  can  be  gained  by  using  12-hourly  heat  fluxes 
that  include  the  effects  of  atmospheric  stability.  This  can  be  achieved  by  using  operational  products 
and  air-sea  temperature  dependent  exchange  coefficients  in  the  NLOM  surface  heat  flux  forcing. 
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Given  that  the  surface  heat  flux  forcing  applied  here  for  the  NLOM  wind  experiments  is  found 
to  cause  the  model  SST  to  relax  towards  an  equilibrium  state  prescribed  by  the  air  temperature 
and  humidity  climatologies,  it  is  reasonable  to  expect  much  better  NLOM  predictions  of  the  SST. 
While  this  is  inadequate  if  one’s  interest  is  simulating  SST  for  climate  studies,  it  serves  well  for  an 
SST  forecast  model,  such  as  the  NLOM,  that  is  used  to  support  Navy  Fleet  operations. 
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Appendix  A 


OXYGEN-BASED  DIAPYCNAL  MIXING  SCHEME 


The  global  mixing  correction  scale  factors  £lk  are  set  to  rates  of  cross  isopycnal  mixing,  based 
on  oxygen  saturation  values  S  calculated  at  mean  layer  interface  depths  zk.  The  rates  of  cross 
isopycnal  mixing  are  determined  in  the  same  manner  as  Shriver  and  Hurlburt  (1997),  and  the 
reader  is  referred  to  the  latter  for  the  physical  rationale.  In  departure  from  this  original  approach, 
a  different  method  is  used  here  to  blend  the  downward  and  upward  vertical  mixing  rates  that 
takes  better  account  of  the  variation  in  oxygen  saturation  with  depth.  Here  the  depth  (Dmj„)  of 
the  absolute  minimum  in  oxygen  saturation  ( Smin )  is  determined  at  each  location,  along  with  the 
greater  depth  ( Dbase )  at  which  the  oxygen  saturation  has  increased  to  5m,n  -I-  1.  These  depths 
provide  a  better  characterization  of  how  the  oxygen  saturation  varies  with  depth.  They  are  used  to 
linearly  weight  with  depth  the  relative  contributions  of  downward  vertical  mixing  from  the  surface 
{wd)-,  and  upward  vertical  mixing  from  the  deep  ocean  (wu). 


The  vertical  mixing  rates  Wd  and  wu  are  the  same  as  in  Shriver  and  Hurlburt  (1997),  and  are 
repeated  here  to  correct  a  couple  of  typographic  errors. 


wd(t)  = 


(e-aDo  _  e-azk 

(1  -  S/So) 


Wu(t)  =  - 


^g-0^2500  _  f>~ CCZk 

(1  -  5/52500 ) 


(Al) 

(A2) 


The  notation  is  the  same,  except  for  the  addition  of  Do  and  So  for  the  surface  depth  and  oxygen 
saturation  value,  respectively.  Here  Do  =  0  m,  D2500  =  2500  m,  and  a  =  3.5  x  10-3m-1  also. 
The  rates  of  vertical  mixing  are  weighted  across  the  region  of  the  absolute  minimum  in  the  oxygen 
saturation  profile  such  that 


Q,k 


Wt d  ^  Dmm 

<  [(1  -  7 k)wd  +  7 kWu]  Dmin  <zk<  Dbase  , 
k  W u  Zk  >  Dbase 


(A3) 


where  the  weighting  factor  is  jk  =  (zk - Dmtn )/(Dbase- Dmm ) .  With  this  weighting  choice,  irregular 
regions  of  high  mixing  rates  are  removed  in  the  fields  obtained  with  the  original  method,  and  the 
fields  vary  smoothly  with  latitude  and  longitude  in  a  manner  consistent  with  the  oxygen  saturation 
climatology  used  (Levitus  and  Boyer  1994).  This  results  in  better  global  detrainment  being  applied 
in  the  NLOM  simulations. 
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Appendix  B 


NLOM  NOTATION 


The  definitions  of  the  symbols  appearing  in  the  many  equations  of  this  report  are  collected 
together  below  for  ease  of  reference. 

Bl.  Layer  Model 

The  quantities  appearing  in  the  equations  for  the  n  active-layer  reduced-gravity  thermodynamic 
NLOM  are  defined  below.  The  units  used  in  the  model  are  indicated  in  parentheses. 

Cp  specific  heat  capacity  for  water,  cf  Table  2 
v*  kth  layer  velocity  (cm  s-1) 

Pk  kth  layer  density  (at) 

hfr  kth  layer  thickness  (m) 

Hk  kth  layer  thickness  at  rest,  cf  Table  3 
Hn  constant  depth  of  ocean  (6500  m) 
pn  constant  rest  density  (27.71  at) 

f  Coriolis  parameter  (s-1) 

po  constant  reference  density,  cf  Table  2 
rw  wind  stress 

Ck  coefficient  of  interfacial  friction 

Ci  coefficient  of  bottom  friction 

u)kt  entrainment  mixing  velocity  due  to  turbulence  (m/s) 

Qk  kth  interface  reference  (ms-1)  vertical  mixing  velocity 
hp  kth  layer  thickness  at  which  entrainment  starts  (m) 

hfc  kth  layer  thickness  at  which  detrainment  starts  (m) 

fife  kth  interface  global  mixing  correction  scale  factor 
Cm  coefficient  additional  interfacial  friction  associated  with  entrainment 
Ah  coefficient  of  horizontal  eddy  viscosity  (m2s-1) 

Kh  coefficient  of  horizontal  density  diffusivity  (m2  s-1) 

Op  reference  coefficient  of  density  climatology  relaxation 
Pk  kth  layer  density  climatology 

Tk  kth  layer  temperature  climatology 

H0  constant  reference  density 

A p~p  minimum  density  contrast  for  (pk+ 1  —  pk) 
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a(T)  coefficient  of  thermal  expansion  for  sea  water  (°C  ]  ) 

r  _  I  g(pn+l  -  Pk)  /  Po  l<k 
kl  \  g{pn+i- Pi)/p0  i>k 


Tk  = 


Wfc  = 


Tw  k-0 

Ckpo\vk  -  Vfc+il  X  (v*  -  vfe+i)  k  =  1, 2, ...  n 
0  k  =  0,  n 

(uifr  +  —  k  1,  2, ...  71 


+12 


u'jf  =  a)fc[max(0,  h £  -  hk)/h 
=  Wfc[max(0,%  ~h^)/h^]2 

_j_  f  0  h\  >  hm 

Um  ~  |  max(0,  Wfet)  hi  <  hm 

Wfc  =  //«f-  uff)dSI  J  J  ClkdS ;  S(z,y)  =  horizontal  area  of  model  domain 
Apfe  =  max(pfc+i  -  Pfc  -  Ap^ ,  0) 


B2.  Mixed  Layer  Model 


The  mixed  layer  equations  of  the  model  use  the  following  notation: 

Tm  temperature  of  mixed  layer  ( °C) 

Tb  temperature  beneath  the  base  of  the  mixed  layer  ( °C) 

hm  mixed  layer  depth  (m) 

Tfc  kth  layer  temperature  (°C) 

Tk(x,y)  kth  layer  temperature  dim.  (°C) 

u)m  vertical  velocity  at  base  of  mixed  layer  (ms-1) 


B3.  Turbulence  Model 

Enclosed  below  is  the  notation  for  the  turbulence  model  equations: 

mi,n  turbulence  model  parameters,  cf  Table  2 
f+  Coriolis  parameter  at  5°  latitude  (s_1) 

h+  minimum  mixed  layer  depth  (m) 

h*m  target  mixed  layer  depth,  >  K 

AT  temperature  contrast  at  base  of  mixed  layer  ( °C) 
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B4.  Surface  Heat  Fluxes 


The  notation  for  the  terms  appearing  in  the  surface  heat  flux  equations  are  given  here. 

Cs  empirical  coefficient  for  sensible  heat  flux 

Cpa  specific  heat  capacity  for  air,  cf  Table  6 

Clo  intercept  value  for  Cl  coefficient 

b  slope  value  for  Cl  coefficient  (°C_1) 

e  ratio  of  molecular  weights  of  water  vapor  to  dry  air 

Mw/Md 

es  saturated  water-vapor  pressure  (mb) 

L  latent  heat  of  evaporation  of  water,  cf  Table  6 

Pa  air  pressure,  cf  Table  6 

@(rr)  Heaviside  step  function 

Ta  climatological  air  temperature  (°C) 

qa  climatological  specific  humidity  of  air 

pa  density  of  air,  cf  Table  6 

scalar  wind  speed  (m/s) 

vmin  minimum  scalar  wind  speed,  cf  Table  6 


Appendix  C 


HEAT  EQUATION 


Heat  energy  within  the  NLOM  is  determined  using  a  layer  averaged  form  of  the  thermal  energy 
equation.  This  equation  is  derived  here  because  it  differs  in  several  respects  from  those  used  by 
other  authors  (McCreary  et  al.  1993). 

Under  the  Boussinesq  approximation,  as  used  in  the  NLOM,  the  thermal  energy  equation  for 
a  fluid  is  given  by  (Kundu  1990) 

^  =  -V-q,  (Cl) 

where  T  is  temperature,  p  is  the  density,  Cp  is  the  specific  heat  capacity  of  sea  water,  q  is  the  heat 
flux,  and  D/Dt  —  d/dt  +  v  •  V  is  the  material  time  derivative.  Viscous  dissipation  of  energy  is 
negligible  in  this  approximation.  We  consider  the  heat  flux  to  be  the  sum  of  a  source  at  the  air-sea 
interface,  and  an  ocean  component  that  obeys  the  Fourier  law, 

V  •  q  =  -Q6{z)  -  V  •  K  •  VT,  (C2) 

Term  1  Term  2  Term  3 


where  Term  1  is  the  total  heat  flux,  Term  2  the  surface  heat  flux,  and  Term  3  the  heat  dissipation. 
The  thermal  conductivity  K,  is  a  second-order  tensor  whose  elements  depend  on  the  local  state 
of  the  fluid.  This  heat  dissipation  term  is  similar  in  form  to  the  viscous  diffusion  of  momentum 
considered  by  Wajsowicz  (1993).  Following  his  same  arguments,  the  tensor  simplifies  under  the 
condition  of  transverse  isotropy  to  a  form  having  separate  horizontal  and  vertical  dependences  (cf 
eq(ll)  in  Wajsowicz  (1993)).  This  yields  after  some  rearrangement 


-PCp  (v  •  VhT  -  +  Q8(z ) 

+V„  •  (K„VeT)  +  ^  (kv|j)  . 


(C3) 


where  VH  is  the  horizontal  derivative,  and  z  is  the  vertical  coordinate.  The  latter  is  defined  to  in¬ 
crease  positively  with  depth  and  is  the  reason  the  vertical  advection  term  has  opposite  sign.  Typical 
values  used  in  ocean  models  for  the  vertical  and  horizontal  thermal  diffusivities  are  Kh  =  300  m2  s-1 
(Hurlburt  et  al.  1996)  and  JCy  =  2x  10~6m2s_1  (Martin  1985),  respectively.  Corresponding  typi¬ 
cal  horizontal  and  vertical  length  scales  are  L  =  10  km  and  d  =  10  m.  This  implies  that  the  vertical 
diffusion  is  negligible  relative  to  the  horizontal  as  K,yL2 /ICnd2  ^7x  10-3.  Only  the  contribution 
from  horizontal  diffusion  is  therefore  retained. 
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The  NLOM  equations  for  ocean  basin  scales  are  derived  by  reducing  the  three-dimensional 
primitive  equations  to  those  for  a  two-dimensional  model  of  a  thin  fluid  on  a  spherical  surface. 
They  are  formally  obtained  by  integrating  in  the  radial  direction  over  the  thickness  of  the  spherical 
shell  while  considering  the  fluid  state  variables  to  be  constant.  Here  we  limit  ourselves  for  simplicity 
to  a  Cartesian  coordinate  system  to  derive  the  general  form  of  the  equations,  as  the  details  on  the 
integration  procedure  for  the  spherical  case  can  be  found  elsewhere  (Moore  and  Wallcraft  1997). 

The  heat  equation  for  layer  k  is  obtained  by  vertically  integrating  the  thermal  energy  equation 
over  the  thickness  h k  =  z^  —  zk  of  the  layer.  Care  must  be  taken  when  exchanging  the  order  of 
derivatives  and  integrals  because  of  the  coordinate  dependence  on  the  integral  limits.  The  relation 
for  differentiation  of  a  definite  integral  proves  useful  for  this  purpose  (Gradshteyn  and  Ryzhik  1980). 

d  C ^ 

—  /  f(x,a)dx  = 

U>OL  J 

+  C“’  (C4) 

J<f>(a)  d(X 

Straightforward  application  of  this  relation  yields  the  layer  averaged  expressions  for  all  but  the 
diffusion  term.  The  latter  is  a  special  case  which  we  postpone  until  further  below.  We  denote 
the  diffusion  for  the  time  being  by  DH-  When  integrating  the  vertical  advection  term  a  physical 
condition  is  applied  when  evaluating  ujT  at  the  layer  interfaces  and  z*  .  There  must  be  no 
increase  in  the  internal  energy  of  the  layer  unless  water  of  a  different  temperature  is  added  to  the 
layer.  This  leads  to  the  layer  heat  equation 

Pk^^-+Pkhkvk-VHTk  = 

Si k~-  +  ^r-  +  Pk  [max(0,  ujk)(Tk+i  -  Tk) 

C/p  C/ p 

-)-  max(0,  -0Jfc_i)(Tfc_i  -  Tk)]  (C5) 

in  which  uk  =  dz^/dt  is  defined  as  the  velocity  of  the  interface  between  layers  k  and  k  +  l. 

The  treatment  of  dispersion  terms  for  layer  models  (D/y)  is  difficult  because  a  physically  consis¬ 
tent  form  cannot  be  directly  derived  from  the  primitive  equation  (Shchepetkin  and  O  Biien  1996). 
The  alternative  is  to  use  a  phenomenological  form  that  best  satisfies  the  constraints  used  in  the  gen¬ 
eral  formulation  of  the  three-dimensional  stratified  flow.  The  form  DH  =  KH^n  •  {hkV HTk)  follows 
naturally  from  the  heat  dissipation  term  in  the  primitive  equation,  with  the  spatial  dependence  of 
the  thermal  conductivity  for  the  layer  being  made  the  same  as  the  layer  thickness  Kk  =  Kh  hk .  It 
satisfies  three  of  the  constraints  used  in  the  general  theory:  1)  in  the  conservation  form  of  the  heat 
equation,  the  dispersion  term  must  have  the  form  of  a  divergence  of  a  tensor;  2)  the  dispersion 
must  be  dissipative,  implying  the  dissipation  on  the  righthand  side  of  the  thermal  energy  equation 
has  the  form  of  divergence  of  a  vector  and  a  negatively  defined  source  term;  and  3)  it  depends  on 
the  symmetric  part  of  a  local  tensor.  In  the  case  of  momentum,  the  equivalent  form  of  dispersion 
as  lateral  friction  has  been  shown  to  give  a  tolerable  solution  (Shchepetkin  and  O  Brien  1996).  The 
similarity  of  the  primitive  equations  leads  one  to  expect  equivalent  results  for  the  thermal  energy. 
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Using  this  choice  for  the  dispersion  term  Dh ,  and  the  definition  of  thermal  diffusivity  k  = 
Kn/{pCp),  one  obtains  the  NLOM  heat  equation  as  follows: 

d(hkTk) 

Pk — g~t - 1-  PkhkVk  ■  v*r*  = 

— |-  npk{y h  *  n))Tk 

L'p 

+  pk  [max(0,  u>k){Tk+i  -  Tk) 

+  max(0,  -wjfc_i)(Tfc_i  -  Tk)] .  (C6) 


The  corresponding  equation  for  the  mixed  layer  is  evident  by  inspection. 

d(hmTm)  ,  rj  m 

Pm.  i  1  *  v  H-*-m  — 

1"  Pm  max(0,  cum)(Tf  —  Tm) 

KJp 

+  Kpmiy H  *  (^mV//))Tm.  (C7) 

The  mixed  layer  density  pm  is  determined  from  the  mixed  layer  temperature  via  the  relation 

Pm  ~  p\  ~  &{Tk)po{Tm  —  Ti).  (C8) 


These  equations  are  noted  to  differ  from  those  of  the  MKM  principally  in  the  form  of  the 
advection  and  diffusion  terms. 
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Appendix  D 


GLOSSARY 


The  acronyms  appearing  throughout  this  report  are  listed  below  for  ease  of  reference. 


ATKE 

Available  turbulent  kinetic  energy 

COADS 

Comprehensive  Ocean- Atmosphere  Data  Set 

ECMWF 

European  Centre  for  Medium-Range  Weather  Forecast 

FNMOC 

Fleet  Numerical  Meteorology  and  Oceanography  Center 

HR 

Hellerman  and  Rosenstein 

IPOPS 

Interim  Primary  Oceanographic  Prediction  System 

MKM 

McCreary  and  Kundu  model 

MLD 

Mixed  layer  depth 

NE 

Northeast 

NLOM 

Naval  Research  Laboratory  Layered  Ocean  Model 

NOGAPS 

Navy  Operational  Global  Atmospheric  Prediction  System 

NRL 

Naval  Research  Laboratory 

OGCM 

Ocean  general  circulation  model 

PSD 

Power  spectral  density 

ss 

Skill  score 

SSH 

Sea  surface  height 

SST 

Sea  surface  temperature 

SW 

Southwest 

TKE 

Turbulent  kinetic  energy 

VAWR 

Vector- Averaging  Wind  Recorder 

WHOI 

Woods  Hole  Oceanographic  Institution 
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